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We consider the energy shift for an atom close to a non-magnetic body with a magnetic moment 
coupled to a quantized magnetic field. The corresponding repulsive Casimir-Polder force is obtained 
for a perfect conductor, a metal, a dielectric medium, with dielectric properties modeled by a Drude 
formula, and a superconductor at zero temperature. The dielectric properties of the superconductor 
is obtained by making use of the Mattis-Bardeen linear response theory and we present some useful 
expressions for the low-frequency conductivity. The quantum dynamics with a given initial state is 
discussed in terms of the well-known Weisskopf-Wigner theory and is compared with corresponding 
results for a electric dipole coupling. The results obtained are compatible with a conventional master 
equation approach. In order to illustrate the dependence on geometry and material properties, 
numerical results are presented for the ground state using a two-level approximation. 
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I. INTRODUCTION 



In recent years [ij, great advances in experimental 
techniques have stimulated an intense theoretical as well 
as experimental activity on the Casimir effect [2i]. In 
particular, the study of the force between an atom and 
a bulk material - the Casimir-Polder (CP) force [1, 01 
- at finite temperature has recently become a popular 
subject of research (see e.g. Refs.fl, H, H Si and refer- 
ences therein). The CP force on e.g. a Bose-Einstein 
condensate is, in fact, not only of theoretical interest 
as the thermal corrections are measurable 

UM- CP 

forces play an important role in a variety of processes in 
physical chemistry, atom optics, and cavity QED. This 
includes e.g. nano-technological applications [uill^, the 
possibility of investigating fundamental forces at the sub- 
micrometer scale [l3| and the effort to miniaturize atom 
chips [3, The CP force can, under some circum- 

stances, turn from an attractive character to a repulsive 
one leading to a quantum levitation phenomena. This 
can be achieved using left-handed meta-materials [Tsl or 
if the system under consideration is immersed in a suit- 
able fluid [IB] ■ As we will argue in the present paper one 
can easily obtain repulsive CP forces if one consider mag- 
netic moment transitions instead of electric dipole ones 
without the need of additional ingredients. Such repul- 
sive forces tend, however, to be very small as compared 
to the conventional attractive electric dipole CP forces. 

Experiments reveal that the proximity of atoms to the 
non-magnetic body may introduce several surface-related 
decoherence effects (see e.g. Refs.flTl [Tsl [l9|). Most im- 
portantly, atoms may be expelled from an atom trap due 
to spin flip transitions induced by Johnson-noise currents 
in the material. If the expelled atoms can not re-enter 
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FIG. 1: Schematic picture of the setup considered in our cal- 
culations. An atom inside a microtrap is located in vacuum 
at a distance z away from a non-magnetic slab with thickness 
h. The slab has an infinite extension in the other directions. 
Vacuum is on both sides of the slab. The slab can e.g. be a 
normal conducting metal or a superconducting metal. 



the trap, it is not clear to what extent the system in gen- 
eral can be described in terms of a conventional master 
equation. 

Speciflcally we will consider a multilevel atom close to a 
non-magnetic slab as illustrated in FiglU where the atom 
may be held close to the slab by e.g. a microscopic atom 
trap [20l | . We focus on magnetic moment couplings to the 
electromagnetic fleld but we will also mention, where ap- 
propriate, relevant effects for couplings to electric dipole 
moments. Furthermore, applying the Heisenberg opera- 
tor equations of motion and the rotating-wave approxi- 
mation lead in general to operator ordering problems and 
great care has to be taken in order to avoid mathemat- 
ical inconsistencies (see e.g. Ref . [2l| ) . The well-known 
Weisskopf-Wigner theory [22] will, however, describe the 
physics in a more clear manner and, as we will see, clari- 
fies the role of various approximative procedures. Numer- 
ical results will be presented by making use of a two-level 
approximation. The two-level approximation for electric 
dipole transitions must be handled with some care due to 



2 



the existence of sum rules [23|. For magnetic transitions 
we will, however, not encounter such issues. 

The paper is organized as follows. In the next section 
we outline the theoretical framework. In order to ob- 
tain the energy shift for ground state atoms it is argued 
that one has to go beyond the conventional rotating- wave 
approximation. In Section IIIII we give explicit expres- 
sions for the appropriate Green's functions in the case 
of a semi-infinite slab and magnetic transitions. Ex- 
plicit expressions for the ground state energy shift are 
then obtained for a zero temperature perfect conductor, 
a metallic slab with dielectric properties described by a 
Drude dispersion relation and extensions of it, a dielec- 
tric medium, and for a superconductor as described by 
a weakly coupled BCS superconductor. In the case of 
a superconducting slab we have to reconsider the low- 
frequency dielectric properties in great detail including 
the presence of impurities. The short and long distance 
behavior of the corresponding CP force, as induced by 
a magnetic moment, are considered in detail. Numerical 
results are compared with the corresponding results for 
electric dipole transitions. By defining a suitable rescaled 
and dimensionless CP force, we can easily compare some 
of our results with similar results for CP forces as in- 
duced by electric dipole transitions. In Section |TV] we 
give some final remarks concerning electric dipole transi- 
tions. Technical details concerning the Green's functions 
used for magnetic transitions, their analytical continu- 
ation as well as various explicit expansion for a perfect 
conductor are presented in two appendices. 



II. GENERAL THEORY 

Let us consider a neutral atom at a fixed position r^^. 
The magnetic moment of the atom interacts with the 
quantized magnetic field via a conventional Zeeman cou- 
pling. The total Hamiltonian is then 

a 

+ J dujhuj{\r,uj) ■{{r,uj) + H' , (1) 

where the effective interaction part is 

H' = l«)(/3|Mo/3-B(r^) . (2) 

a 

Here f (r, uj) is an annihilation operator for the quantized 
magnetic field, \a) denotes the atomic state and Ea is the 
corresponding energy. We assume non-degenerate states, 
i.e. Ea Ep for a ^ (3. The magnetic moment of the 
atom is iJbap = (a|/i|/3), where fx is the magnetic moment 
operator. The magnetic field B(r) = B(+'(r) -I- B(~)(r) 
is written B(+)(r) = V x A(+)(r), where B(-)(r) = 



(B(+)(r))^ and where the vector potential is 



X G(r,r',c.')-f(r',c.') . (3) 

Here the imaginary part of the complex permittivity is 
e2{r, uj). The dyadic Green's tensor G(r,r',tt') is the 
unique solution to the Helmholtz equation. Because the 
Helmholtz equation is linear, the associated Green's ten- 
sor can be written as a sum according to 



G(r, r', uj) = G°(r, r', uj) + G'5(r, r', uj) 



(4) 



where G'^{r,r',uj) represents the contribution of the di- 
rect waves from the radiation sources in an unbounded 
medium, which is vacuum in our case, and G^{r,r',uj) 
describes the scattering contribution of multiple reflec- 
tion waves from the body under consideration. The pres- 
ence of the vacuum part G'^{r, r', uj) in Eq.([4]) will in gen- 
eral give rise to divergences in the energy shifts to be cal- 
culated below. A renormalization prescription is there- 
fore required. We subtract the vacuum part G*'(r, r', cj), 
i.e. we neglect possible finite corrections due to this vac- 
uum subtraction like conventional vacuum-induced Lamb 
shifts. Since we, in the end, are going to consider CP 
forces all additive coordinate independent corrections to 
energy shifts will, anyway, not contribute. When we be- 
low refer to a renormalization prescription the procedure 
above is what we then have in mind. 

We now consider the Hamiltonian in Eq. ([T]) and apply 
the well-known Weisskopf-Wigner theory for the transi- 
tions a ^ (3, where Ea > Ep. The solution to the time- 
dependent Schrodinger equation in the rotating-wave ap- 
proximation (RWA), i.e. applying H' w HjiwA, where 



-ffwA = - ^ |a)(/3| Ma/3 • B(+)(rA) -f h.i 



(5) 



I3<a 



is then {uja = Ea/h) 

\m) = Ca{t)e-'^-*\a)^\0) 

d\ d^Y^T. 

•'^ m=lp<a 

X c„0(r,c^,t)e-'("+'^'')*|/3)® |l™(r,c^)) . (6) 

Here the initial state of the atom-field system is |a) (8) |0), 
where |0) denotes the vacuum of the electromagnetic field 
and |lm(r,a;)) = fl^{r,uj) |0) is a one photon state. We 
also make use of the notation /3 < a in a /3-sum to denote 
the condition Ep < Ea- The coefficient Ca{t) is then 
determined by 

dca{t) 



= f dt' K%{t-t')ca{t') , (7) 

•^0 0<a 

where K^p{t) is the renormalized version of the kernel 
1 /"^ 

Kcpit) ^ da;e-^(— °'')*r„^(rA,a;). (8) 



Here we have defined = {Ea — Ej3)/Ti > as well as 

(9) 



^ • Im[^ X G(r, r, fip^ , 

which is the spontaneous spin flip rate for the transition 
a ^ (3 where Ea > Ep. Assuming that the Markov ap- 
proximation holds, i.e. memory effects can be discarded, 
we can make the substitution Ca{t') — > Ca(t) in Eq.(l7]). 
By extending the remaining time integral to infinite time 
and making use of the distributional identity 



the solution to Eq.([7]) is then of the form 

Ca{t) = Cq(0) expi - [ ir„(r^,Wa) 



(10) 



(11) 



The spontaneous spin flip rate for an atom initially in 
the state la) is 



a. 



a. 
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FIG. 2: Feynman diagrams. From left to right: the virtual 
process in vacuum, stimulated emission in the heath bath, and 
the absorption-emission process in the heath bath, all leading 
to the same final state a. 



where 5uja\v,uja) = Y.pKa^^a'pi'^^^aP) = l'^)- Here 
we have defined the resonant contribution 



• Re['^ X G(r,r,w) X ^] • , (16) 



Mo 



ra/3(l",^a/3) 



(12) 



where r^^(r, (jj^^) is the renormalizcd version of the 
spontaneous spin flip rate in Eq.®. In Eq. pdj) we have 
also defined the frequency shift 6uJa{r,LUa) by 



0<a 



(13) 



where again Suj^p{r,uj) is the renormalized version of 

. ^ 1 duj' . M . ^ 

6iJc.0{r,uj) ^ tt'P 7 r^Mr,^') , (14) 

Z TT Jq LO — Ul 

where Tap{r,uj) is given in Eq.Q. Sinnlar results to 
Eqs.® and lU]) are also derived in Ref.'24]. These lat- 
ter results can, of course, be obtained using conventional 
perturbation theory (see e.g. Ref.[25j). 

In passing we observe that Eq.®, i.e. the spontaneous 
spin flip rate, has been derived using the Heisenberg 
equations of motion in Ref . [29| . Despite the fact that 
such an approach leads to a well-known mathematical 
inconsistencies due to operator ordering problems [2l| , it 
gives, nevertheless, the same result. Similar results for 
the frequency shift for the more well-known case of an 
electric dipole can e.g. be found in Refs.[27l. [2^. 

The principal integral in Eq. (jl4[) can be further evalu- 
ated by means of well-known contour integral techniques 
due to the absence of poles in the upper complex fre- 
quency plane (see e.g. Refs. [2l1. [27l [29j 1 . For the transi- 
tions a ^ /3, where Ea > Ep, we obtain 



and the off-resonant contribution 

,(2) 



/•OO 

n- Jo 



^'^a/3(l",t^) = 

, AtQ/3 • \^ X G{v,v,iuj') X ^] • ^Ipa 



(17) 



In the RWA, the energy shift for the ground state in 
the Weisskopf-Wigner approach is, however, vanishing in 
the large time limit. On the other hand, applying the 
interaction H' « Hawr, where 

--^|a)(/3|M„^-B(-)(rA)+h.c. , (18) 

rather than Hjuva when considering the transitions /3 — s- 
a, where Ef^ < Ea, i.e. ujfja < 0, the decay rate vanishes 
but the energy shift is non-zero. It is given by 



Sujf3{rA,UJp) = Sujj^ai^A,i^f3a) , 

a>f3 



(19) 



where again Sujp'^{r,uj) is the renormalized version 
Eg. p^ . In the same fashion as above one may show 
that for the state |/3), the frequency shift is 



(2) 

6cji3{rA,uj0) = 5uj^p '{rA,ojp) , 



(20) 



/2') ^ {'2') 

where now (5a;^ '^(r,W;3) = ^^^^ (5wL (r, w^^). For states 
such that Lopa < there will be only off-resonant contri- 
butions when converting the cj-integration in Eq. (|14p to 
imaginary frequencies. 

From these discussions we now realize that, in general, 
it is not sufficient to make us of only the effective inter- 
action HfiwA in the rotating-wave approximation. In- 
stead we must use the complete interaction Hamiltonian 
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H' as given by Eq.(l2])- This has the unfortunate conse- 
quence that an exact analysis is not any more possible. 
The Weisskopf-Wigner theory can, however, easily be ex- 
tended to a more general situation, at least when one 
limit oneself to self-energy corrections to at most second- 
order in perturbation theory (see e.g. Ref.llBl). The 
energy shift for a given state of interest to us is then ob- 
tained by calculating the second-order energy shift using 
the complete interaction Hamiltonian Eq. If we limit 
ourselves to perturbation theory an alternative procedure 
is to consider the coordinate-dependent energy shift and 
making use of linear response theory in a standard man- 
ner (see e.g. Ref.^3Q]). We are then interested in finding 
an expression for {'il}{t)\H'\tl;{t)) = i{tp{t)\H'j{t)\'ip{t)) i , 
where in the interaction picture and to lowest order in 
H' 



\m)i = mo)) + T- / dt' H',{t')mo)) . (21) 

«^ Jto 

Here we eventually consider the limit Iq — > — oo for the 
initial time ^o- In terms of a retarded response function 
Vjiitjt') we can therefore write 

1 f°° 

m)\H'm))^^ dt'VnM, (22) 



where 



iVH{t,t')^e{t-t'){i;ito)\[H'j{t),H'j{t')]mo)) , (23) 

and where we have assumed that (■!/'(<o)|^f'|V'(*o)) = 0. 
The retarded Green's function Eq. ([23|) is then directly 
related to the retarded Green's functions for our basic 
quantum fields A*^+^(r), and its hermitian conjugated 
field, as given by Eq.Q, and can be calculated in a stan- 
dard manner with a knowledge about the initial state 
lipito)). The virtue of this formulation is that we now can 
allow for any quantum state of the radiation field, like a 
thermal background, of photons by making use of a suit- 
able random phases for a pure state of the radiation field. 
One is then led to a consideration of so called real-time 
finite temperature propagators which has been consid- 
ered elsewhere in great detail (see e.g. Refs.[3l|). In the 
present paper we will, however, restrict us to an initial 
vacuum state of the radiation field and return to the case 
of a background of thermal photons elsewhere Ref . [33| . 
The energy shift then obtained will contain processes cor- 
responding to photon absorption and emission, as well as 
stimulated emission, with a thermal background of pho- 
tons. This physical picture is illustrated by means of the 
Feynman diagrams in Figl^] The initial state of the atom 
is then assumed to be fixed which has been referred to as 
a constrained free energy of an atom immersed in heat 
bath Ref. [13. 

For reasons of completeness we now write down the 
energy shift for a background of thermal photons, at 
a temperature T, including finite temperature sponta- 



neous emission and absorption processes, valid to second- 
order in perturbation theory using the Feynman dia- 
grams as indicated in FigO We obtain (5a;Q(r, Wq, T) = 
J2f3 <^^Q/3(i"i ^a0, T), where (5w^^(r, uj, T) is the renormal- 
ized version of 



1 

— V 

27r 







5ijJal3{Y,UJ,T) = 

- {l + n{u'))Vc,i3{v,J) + 



LU ~ LU' 



27r 



(24) 



where Tajj(y^Lo) again is as given in Eq.([9]) and n(a;) = 
\ / {expihuj / ksT) — 1) is the Planck black-body distribu- 
tion. As the temperature T ^ in Eq. ([M)) we obviously 
have a smooth limit (see in this context e.g. Ref. [8]). 
Similar results to Eqs. ([9]) and lfT4|) are also derived in 
Ref. [tI. [2^. [3^. [35! for electric dipole transitions. 



III. ATOM NEAR A NON-MAGNETIC SLAB 



Let us now consider the geometry as shown in Fig. [U 
i.e. a semi-infinite slab with finite thickness. The ana- 
lytical continuation of the equal position Green's tensor 
Eq. (|H|l is then 



X G^{Y,Y,iuj) X = 



1 

47r 



/|| (cj) 
/ii(^) 
I±[uj) 



(25) 



where to is real and 
1 



^11 H 



X 







?7o(A,a;) 



- 2 770 (A.cj) js 



OJ 



Cjv(A,tj) -I- r/o(A,w)Cji/(A,w) 



(26) 



'7o(A,w) 



'2"«(^^'^)^CM(A,tj). (27) 



The scattering coefficients Cn{X,u!) and Cm(A, w), which 
are analytical continuations of Eqs. (|A2p and (|A3[) . are 
given by [36|] 



CAr(A, w) = rp{X,uj) 
Cm{X,uj) = rs{X,uj) 



- 2ri{\,Lu) h 



l-r2(A,w) e-2';(>>'^)'« ' 

2^ _ 2 J)(A,w) h 

l-r2(A,w) e-2'?(A,^)ft ' 



(28) 
(29) 



with the electromagnetic field polarization dependent 
Fresnel coefficients 



rs{X,uj) = 
rp{X,uj) = 



?7o(A,cj) ~ i]{X,uj) 
r7o(A,cj) +?7(A,w) ' 
e{iuj)i]o{X,uj) - r]{X,oj) 
e{iuj) 77o(A, Lj) + 77(A, lu) ' 



(30) 
(31) 
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which are analytical continuations of Eqs 
(|A5p . Here we have defined ri{X,uj 
770 (A, = V/c^ + and k 



1) and 



v/^MH + A^, 
Lo/c. Below we will mostly 



consider the large h limit, i.e. h is supposed to be large 
in comparison with any other length scale in the system. 

Due to causality, the complex dielectric function 
e{oj) = ei(w) + 162 (w) will in general obey the Kramers- 
Kronig relations. Hence [s^l 



e[tuj) = 1 H / ax ■ 



(32) 



where e^iuj) is a real function. Due to physically neces- 
sary conditions, 62(0;) > for a; > (see e.g. Ref.[37|). 
This guarantees, in general, that G'^(r, r', iu) is real. 

The Casimir-Polder force obtained from any frequency 
shift in 6LUfi{r, up) is now given by 



(33) 



which actually is independent of the additive and di- 
vergent contribution due to the free Green's function 
G°(r,r,o.). 

Let us also limit our attention to ground state transi- 
tions /3 — > a, where Ep < Ea or copa ^ 0, in which case 
we may apply Eqs. ([20|) and (|25|) . Only the off-resonant 
term contribute in this case, i.e. 

Su;p{z,u;p) = £| ^ j(|^f"|2 + ^^^) 

+ |Mfpzj_(z,c.0„) 1 (34) 



(35) 



where we have defined ( using p — 1 1 , _L ) 
u dw_ 



TT 







In the numerical analysis we have found it useful to 
make use of the following property of ip(z,uj), i.e. we 
can write 



ip(z,w) = - 



dx 



l + x 



Ip{ujx) 



2 -'P 



dx 



Ipiuj/x) , 



(36) 



with a natural choice e = 1. As a curiosity we remark 
that such a behavior of integrals under a modular trans- 
formation X 1/x has been noticed before and found to 
be useful in e.g. quantum electrodynamics (38| . 

A. Perfect conductor 



The integrals Eq. (|35|) can, in general, not be computed 
analytically. For a perfect conductor Cat (A) = 1 and 
Ca/ (A) = —1 for any h. The corresponding integrals can 
then be simplified. In this case, the integrals Eqs. 
and ([27)1 are related by 



1 



(37) 



For an atom in the state where Ep < Ea, it 
is straight forward to show that the frequency shift is 

(kpa = \(^Pa/c\ = kap): 



inix) 
ij_{x) 



+ l^f pz±(fcc./3^) 1, (38) 



(39) 
(40) 



5ujp{z,ujp) 



where the dimensionless integrals are given by 



d^ 



2x 

ix 



IT J, i2xy+e 



This result is analogous to Eq.(10.87) in Ref.[23, which 
describes the attractive potential for an electric dipole 
outside a perfectly conducting plate. The magnetic fre- 
quency shift for the excited state is given in the Appendix 
[M see Ea. (IBl|) . We observe that the energy shift corre- 
sponding to Eq. (j38p is independent of h, i.e. it may be 
derived purely in classical manner. 

For short distances (non-retarded limit), i.e. kapz <^ 
1, one may easily show that i±{kapz) « 2i^\(kapz) w 1, 
in which case Eg. ((38)) is reduced to 

Mo 1 



Sujp{z,LL!p) 



647r h z^ 



E 



Mo 



647r hz^ 



\fi- fi + p.^p^\ 



(41) 



due to completeness of states. Numerical studies show 
that Ea. (|¥T)) actually a good approximation for kapz < 
0.01. This result may also be obtained using a classi- 
cal approach, e.g. the method of images for a mag- 
netic moment (see e.g. Ref . [39l| ) . The frequency shift 
in Eq. ([1T|) corresponding to the repulsive Casimir-Polder 
force Fp{z,ujp) as given by Eg.pS]) i.e. 



Fpiz,u;p) 



3^0 
GAttz^ 



(42) 



Here is a unit vector in the z-direction, i.e. normal 
direction to the plane of the conductor. With a = x,y,z 
and fi'^f^ = gseUS^^ /2me, where S^^ = {a\Sa\P) are 
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dimensionless, the frequency shift ^(^^(zjW/Jq) as given 
by Eq. ([lT|) can be re- written as 



9l 



647r ei^hz 



\S-S + S,S,\P), (43) 



where Ae = fi/nieC is the Compton wavelength. For 
electric dipole transitions, the corresponding frequency 
shift can be written in an analogous form where the 
magnetic dipole moment /Xq/j is replaced by the elec- 
tric dipole moment d^^ = er^^ (see e.g. Ref.|27|'). 
where Yap = (ajfl/?). Here r is the displacement vec- 
tor, whose order of magnitude is the Bohr radius aq. As 
(Ae/ao)^ ~ 10^^ we realize that magnetic dipole transi- 
tion CP forces are in general weak as compared to the 
corresponding electric dipole transitions CP forces. 

In the long distance (i.e. retarded) limit, correspond- 
ing to kapz ^ 1, we realize that i±{ka(3z) ~ i\\{ka(3z) ~ 
2/{tt kapz). Numerical studies show that this is a good 
approximation for kapz ^ 10. Eq. (|38ll is then reduced to 



6ujf3{z,LL!i3) 



Mo 



167r2 hz'^ 



y — 



,/3a|2 



(44) 



and the corresponding repulsive CP force will then have 
a 1/z^ dependence. Eq. (P^ is an complete analogy with 
the famous Casimir-Polder energy shift for electric dipole 
transitions in the large distance limit 0]. The factor 
containing the sum over intermediate states in Eq. (l44p 
is then related to the static polarizability of the atom 
considered. In our case the corresponding factor could 
similarly be interpreted as a static magnetization of the 
atom. 



Metal 



In vacuum, the Drude dispersion relation may be writ- 
ten (see e.g. Ref. 40,]) 



e{iuj) = 1 + 



UJ (uj + u) ' 



(45) 



where for gold the plasma frequency is cUp = 9.0 eV and 
the relaxation frequency is = 35 mcV. 

Let us first consider the case v — and hence e(icij) = 
1 -|- LUp/uj'^. In order to present analytical as well as nu- 
merical results we now consider a two-level approxima- 
tion {P = g and e = a) for an infinitely thick slab. The 
magnitude of the ground state CP force is 



Fg{z,ujA) = -dSujg{z,ujA)/dz , 



(46) 



with uja = ^e~^g- It is now convenient to explicitly per- 
form the spatial z— derivative in equations Eq. (j26p and 



(|27|) . This enables us to define a rescaled and dimension- 
less force quantity Fm{z.,uja) = FAiikAz) for magnetic 
moment transitions as 



FM{kAz) = Fg{z,u!A)- 



■yi-Kz" 



= {\S^\^ ^\Sy\^)\{kAz) + \S,\'-~%y_{kAz) , (47) 

where Sa = , a = x,y,x, and fc^ = uj/c. Here we 
have defined the functions 



ipikAz) = 



^ Jo {kAZ?+e ' 
where p = || or _L and a = ujp/ujA^ in terms of 

J_,(^,^) = 24 Hdxx'e-'^^ 
Jo 



IpitakAz) , (48) 



(49) 



as well as 



/|l (^, 7) = 2M dxxe 



y.\{e+x') 



v/t^TcM^ -I- VC^Tx^ 



{e + i^we cVt' +e+x^ 



(50) 



Even though the Eqs. (l49l) and ((50|) appear complicated 
their asymptotic expansions can, however, be obtained in 
a straightforward manner. If fc^z ^ 1 then, according to 
Eg. (1^5]) . only small ^ contributes. If, in addition, akAZ ^ 
1 then 



I±iC,akAz) ~ {akAzY 
J||(^, afcAz) ~ (afc^z)^! i + 



2e 



2 (akAzY + 2e 



, (51) 



where we refrain from writing down the higher order 
terms. 

If, on the other hand, akAZ ^ 1 we find that 
Ip{£,, akAz), for a = ||, _L, does not depend on a, and 

/_L (e, akAz) ~ 2(3 6^ -f 4^2)6-2? ^ (53) 

as well as 

/,, (5, akAz) ~ (3 6^ + + 8C^)e-2« . 



(53) 



Eqs. ([47| and (|48|) then imply that if kAZ and akAZ ^ 1 
then 



1 



FM{kAz)^\{\S^\' + \Sy\')\- + 



2 2 + \/2o 
,2 \ {akAzf 



(54) 
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IQ-" L 1 i_iJ 1 1^ ; i_J ; ^1 ; ^1 i 

10-^ 10"^ 10 2 10-1 J j^g2 ^q3 

kAZ 

FIG. 3: The ground state dimensionless CP force FM{kAz), 
as defined in the main text Eq. (|47[) . in the case of a dielectric 
constant e(iLj) = l + w^a'^/tj^ with a = lop/loa, as a function 
of feA2 = u>Az/c. For convenience we have chosen \Sx\^ = 
1 5*5,1^ = \Sz\^ ~ 1/4. The curves correspond to a = 1 (lower 
solid curve), a = 10, 10^, 10^, 10"*, and a = oo (upper curve). 



10"'' 





, 


/// y y ' 
-/// //'' y' eiiiS) 


= 1 + ( ± ^ - 


'O-'v = 10^ 





10-* 10-^ 10-2 ^ ;^q2 



FIG. 4: The ground state dimensionless CP force FM(kAz), 
as defined in the main text Eq. (|47p . in the case of a dielectric 
constant e{iuj) = \+ijj'^/ijj{ll} + v) with uja = i^p, as a function 
of kAZ — LOAzjc. For convenience we have chosen \Sj^ = 
ISyp — \Sz^ = 1/4. With V = v/loa the curves correspond 
to V — Q (upper solid curve), v = 1, 10, 10^, 10^, and v = lO'' 
(lower curve). 



i.e. Fg{z,u)A) oc a'^/{kAzY- If, on the other hand, 
akAZ <^ 1 but kAZ sufficiently small, then 

FM{kAz) ~ ^ + \Sy\^ + 2|5,|2) , (55) 

i.e. Fg{z,ujA) oc l/(fcyiz)^. With akAZ ^ 1 and fc^z > 
10, we obtain 

FMikAz) ~ (I^.P + \Sy\' + \S,\') , (56) 

ttkaz 

i.e. Fg{z,ujA) oc l/(fc^z)^. In FigO the scaled and di- 
mensionless CP force F^iikAz) is presented and the nu- 
merical calculations are in excellent agreement with the 
asymptotic expansions as given above. The choice for 
the matrix elements \Sx\^ = |S'yp = \Sz\^ = 1/4 is not 
important and we do not in general find any essential 
difference with regards to r- and p-polarizations. We will 
therefore, for convenience, make use of the same set of 
matrix elements in the numerical calculations throughout 
the present paper. 

Now returning to the Drudc formula Eq. ((45)) with a fi- 
nite ly, we observe that the a = Wp/cj^-parameter in the 
Eqs.lgg]) and ((501) will be rescaled, i.e. + i^) 

with 9 = v/loa- For the parameters for gold mentioned 
above and for rubidium atoms with lja = 27r x 560 kHz 
we then find that a « 3.9 • 10^ and v k, l.h ■ 10^ and 
one therefore expect that we are close the perfect con- 
ductor limit. A numerical investigation verifies indeed 
that this is the case and Fj^ikAz) is practically indistin- 
guishable from the perfect conductor limit. For reasons 
of completeness we show FM{kAz) in FigH] for a = 1 
for various values of v. It is then clear that it is only 



for a not to large that we are able to see the effect of 
dissipation due to the presence of a finite v. 

More elaborate experimental investigations of the op- 
tical properties of gold have been carried out by Decca et 
al. 41] and by Svetovoy et al. [i^]. In Ref.filj it was, in 
particular, argued that measured Casimir pressures could 
be used to find constraints on possible Yukawa couplings 




kAZ 

FIG. 5: The ground state dimensionless CP force FM{kAz), 
as defined in the main text, in the case of a dielectric constant 
e(iLj) — l + LOp/uJ^ as a function of kAZ = ujaz/c (lower curve). 
For convenience we have chosen \Sx'^ ~ \Sy'^ = \Sz'^ = 1/4. 
The upper curve shows FMikAz) when including interband 
core electron transitions according to Ref. |4l|] in terms of 
Ea.(|57p with parameters as quoted in the main text. 



8 




FIG. 6: The ground state dimensionless CP force FmikAz), 
corresponding to Eq. (|34|l in the main text, using a two-level 
approximation. We make use of a dielectric constant tiiuj) — 
1 + io^a^ / u)'^ , with a = ujp/uJA, as a function of kAZ = loaz/c 
for a fixed value a = 10* . For convenience we have chosen 
\Sxf = \Syf = \Sz\'^ ^ 1/4. The curves correspond to h/z = 
oo (upper solid curve), h/z = 10^^, 10~^, 10"*, and h/z = 
10"'' (lower curve). For h/z 
the infinite slab limit. 



0(1) we are already close to 



from non-standard particle physics. In this work a gener- 
alized 6-oscillator form of the dielectric, constant taking 
into account interband transitions of core electrons, was 
given in the form 



e(iw) 



6 

E 



(57) 



with Up = 8.9 eV and where the rest of the parameters 
fj and gj were fitted to known tabulated optical data. 
We have used the corresponding set of parameters in a 
calculation of rescaled and dimensionless force quantity 
FuikAz) in Eg. ([17)) in order to reveal the effect of a 
more precise description of the dielectric properties of 
gold. The results are presented in FigOand we conclude 
again that in general it will only be for loa — Wp, that 
details in the electromagnetic response properties of gold 
will be revealed in FM{kAz) as defined by Ea. (|47)) . 

In Fig[6]we illustrate the dependence of the thickness 
h of the slab on ground state dimensionless CP force 
FuikAz) in the two-level approximation, corresponding 
to Ea. ([M)) . for an almost perfect conductor [up/uJA ^ !)■ 
We observe that for h/z = 0(1) we are already close to 
the situation of a infinitely thick slab. We also notice 
that the dependence on the geometrical ratio h/z is very 
similar to the dependence on the material property ratio 



LOp/ujA for an infinitely thick slab as in FigO 
C. Dielectric medium 



As an example of a dielectric media we consider sap- 
phire for which the electric dipole induced CP forces has 
been considered in Ref.ji^. We find that the typical 
"two-plateau" behavior of e{iu}) of sapphire as obtained 
numerically in Ref.[43| can be parameterized in the fol- 
lowing form 



e{iuj) = 1 + 



,2 ' 



(58) 



with the parameters iHpi — 0.16 eV, L0p2 = 30.8 eV, wi = 
0.07 eV, L02 = 20.8 eV, where lu is given in units of eV. In 
FiglT] we illustrate the behavior of the ground state two- 
level rescaled force FM{kAz) where we now observe that 
Fg{z,u:A) oc l/{kAzY even in the case rubidium atoms 
with a very small uja = 27rx 560 kHz as compared to other 
parameters in Eg. ([55]) . This change of slope in FM{kAz) 
as a function of kAZ as compared to a metal is due to 
the fact that all the parameters in Eg. fSS]) are large as 
compared to oja- Indeed, using the dielectric constant 
Eq.dHHl) and Eq.dJT]) as well as Eq.fgS]), one finds that if 
kAZ <^ 1 then ip{kAz) ~ kAZ for p = || or _L. 




FIG. 7: The ground state dimensionless CP force FMikAz), 
as defined in the main text, in the case of sapphire with a two- 
plateau dielectric function as given in the main text Ea. (|58ll . 
For convenience we have chosen \Sx^ ~ \Sy^ = \Sz'^ = 1/4. 
The upper dashed curve shows FuikAz) when uoa = ckA = 
2ti X 560 kHz. The middle dotted curve corresponds to lua = 
uop for sapphire. For comparison we also give FM(kAz) with 
a dielectric constant e{iLjj) = 1 + ijj^/oJ^ with uja = (lower 
solid curve). 
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D. Superconductor 



as well as 



As the total current density is assumed to respond lin- 
early and locally to the electric field, the dielectric func- 
tion e{Lu) at zero temperature can be written in the form 



e[uj) = 1 h i 



eouj 



(59) 



Here cr(w) = cri(cL;) -I- i(72i^) is the, in general fre- 
quency dependent, complex optical conductivity. A de- 
tailed and often used description of the electrodynamic 
properties of superconductors was developed by Mattis- 
Bardeen and independently by Abrikosov-Gor'kov- 
Khalatnikov [i^, based on the weak-coupling BCS theory 
of superconductors. In the clean limit, i.e. I ^ ^Qj where 
I is the electron mean free path and is the coherence 
length of a pure material, the complex conductivity, nor- 
malized to cr„ = (Ti(T — Tc), can be expressed in the 
convenient form 14611 



-huj 



where 



Ui{x)u2{x) 

Here we have defined the functions 



(60) 



(61) 



U2{x) = ^J{x + hiuy 



A2 



(63) 



Here, the well-known BCS superconducting energy gap 
A is given by A = 3.53 kBTc/2 in terms of a critical 
temperature Tc [s^- In Fig[5]we plot ai{uj) and 172(0;) as 
functions of the dimensionless variable q = A/Tiui. Since 
the low-frequency properties of e^iu) will dominate the 
contribution to CP forces, it is important to notice that 
Eg. ([60)1 leads to the asymptotic behavior (T2(a;)/(J„ 
ttA/Tiuj as huj/A 0. The dielectric function e((jj) will 
therefore have a pole of second-order at w = 0. 

The Kramers- Kronig relation Eq.((32| has then to be 
modified according to a standard procedure [37j with the 
result 



e{tco) = l+^ + -P 



dx 



X £2(2:) 



(64) 



where we have defined = 7r(T„A/eo?i- For niobium 
hujsp = 2.4 eV which can be compared to TicUp = 9.0 eV 
for gold. 

As shown by Anderson and others [13,14^, the presence 
of non-magnetic impurities, which we only consider in 
the present paper, will not modify the superconducting 
energy gap as given by A. The complex conductivity will, 
however, in general be modified due to the presence of 
such impurities. In terms of the normal conducting state 
relaxation time r, due to the presence of non-magnetic 
impurities, one can write 





0-2(9) / 






cri(g) / 













0.5 



1 



1.5 



FIG. 8: The complex conductivity according to the Mattis- 
Bardeen theory. The variable q is given by g = A/hw. We 
also observe that ai(uj) = if g > 1/2. The solid curves 
correspond to an infinite relaxation time r. The dotted curves 
correspond to niobium for which we choose 13.61 = R/tA, 
which corresponds to a mean free path of 9 nm. 



sgn{x)^/x'^ - A2 if |A| < a; , 

ui{x) = { (62) 

-iVA2 -a;2 if|A|>a;, 



0-2(0;) -I- i(Ji{uj) 



1 

2UJT 

1 

~2WT 



A 



dx 



lix) + 1 



gix) - 1 



dx 



A 



U2 ~ ui + ih/r 
9{x) - 1 



U2 -|- ui — iTi/t 

g{x) - 1 



U2 + Ml — iTi/t U2 + Ul + iTi/t 



(65) 



For niobium we choose r such that Ti/tA — tt^q/I — 
13.61, corresponding to the experimental coherence 
length ^0 — 39 nm and the mean free path Z(r « 9K) = 
vpT = 9 nm [50]. The low- frequency limit of Ea. (|65p can 
be calculated in a straightforward manner and we find 
that 



cr„ Tllu \ At 



(66) 



as huj/A 
function 



0, where we have defined the real- valued 



fix) 



In (a;/2 + ^{x/2y - 1 ) 



1 



(67) 



We notice that in the presence of non-magnetic impuri- 
ties Eq. (|M|) we still obtain a l/w singularity in cr(w) for 
small uj. In Figl5]we illustrate the w-dependence of (t{uj) 
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in the clean limit (r — > cx)) as well as for niobium. Nu- 
merical investigations show that the asymptotic formula 
Eg. ([66]) actually describes cj2{ijj) for a wide range of fre- 
quencies. We also notice that for any finite t/Ti/S. ^ we 
have the asymptotic behavior ct(w) — > as Tilj/ IS. — > oo. 
The discussion in this Section reveals that for supercon- 
ductors, we are led back to the considerations for a metal 
as discussed in Section IIII Bl apart from the continuum 
contribution in Ea. ([64)l . A closer of investigation shows 
that its only for frequencies uua ^ that the CP force 
will be sensitive to this continuum contribution. For such 
frequencies, however, the Cooper pairs of the supercon- 
ductor start to break up and we are again led back to a 
metal. 




10-" 10-^ 10-2 jQ-i ^ JO 10^ 10'' 



FIG. 9: The attractive ground state dimensionless CP force 
FsikAz), as defined in the main text Eq. (|68p for electric 
dipole transitions. For convenience we have chosen \dx\'^ — 
\dy\^ = |dz|^ = 1/4. With a dielectric constant eiiuj) — 
1 -l-oip/w^, the middle curve shows \FE{kAz)\ as a function of 
kAZ = uoaz/c with to a = i-Op and, for comparison, the lower 
curve corresponds to FMikAz) (repulsive) for uja ~ i^p and 
IS^P = \Sy\'^ = \Szf = 1/4. The upper curve shows the 
same FM{kAz) as well \FE{kAz)\ in the perfect conductor 
limit uip ^ oo. 



IV. FINAL REMARKS 

In addition to the interaction Hamiltonian H' = — /x-B 
as given Eg . Q we may also a consider the Power- Zienau- 
Wooley [5l|, \^ Hamiltonian i/' = — d • E in order to 
describe electric dipole transitions. For the semi-infinite 
half space geometry as considered in Section IIIII one can 
now easily transfer the results for magnetic transitions 
to electric dipole transitions according to the observa- 
tions in Appendix |^ This has already been alluded to 
in Section IIII Al There is, however, a subtle difference 
due to the role of the scattering coefficients Cat (A, us) and 
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10-" 10-^ 10-2 jQ-i 1 10 10^ 10^ 



kAZ 

FIG. 10: The absolute value \FE{kAz)\ of the attractive 
ground state dimensionless CP force FsikAz), as defined 
in the main text Ea. (|68|l for electric dipole transitions with 
\dx^ ~ \dy^ — |rfzp = 1/4. The upper dashed curve gives 
\FE{kAz)\ for the rubidium transition lua = ckA = 27r x 560 
kHz and in the case of sapphire with the two-plateau dielectric 
constant as given in the main text Eg. dSSp . For comparison 
the middle solid curve shows \FE(kAz)\ for a dielectric con- 
stant e{iuj) = 1-1- Lo^/uP with to a = ujp. The lower dotted 
curve shows \FE{kAz)\ in the case of sapphire but now with 
LOA = 9.0 eV. 

Cm(A,u;) in Egs. ([26|) and ([27| which, for electric dipole 
interactions, has to be interchanged. We then define a 
dimensionless and rescaled attractive ground state CP 
force FEikAz) in the two-level approximation analogous 
to Eq.dlTl), i.e. 

= + |4nf (fc^z) + lll^lfikAz) , 

(68) 

where F^{z,uja) is the ground state CP force induced 

by electric dipole transitions and \da\'^ = |daP/|dp, a = 
x,y,z. The functions if_ and correspond to Eq. ((48|) 
with the only difference that the scattering coefficients 
CAr(A, cj) and Cm(A,cl;) have been interchanged. Similar 
to the analysis in Section IIIIBl i.e. for a metal, and for 
kAZ and akAZ ^ 1, one may show that 

f(fc^z)c.-|— ^ , 7f(fc^z)^--|^ . (69) 

If, on the other hand, afc^z ^ 1 and fc^z ^ 10 we obtain, 
apart from a minus sign, the same asymptotic expansion 
as in Eq. ([5S|) by replacing jS'^p with \da!^ for a — x,y, z. 
In Figl9]we illustrate these limits for FE{kAz) where we 
observe that it is only for magnetic moment transitions 
we obtain a behavior Fg(z,ujA) oc /(kAz)'^ for small 
sufficiently kAZ. 
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In Fig [To] we show the attractive CP force FE{kAz) 
for sapphire as well as for a metal and for different fre- 
quencies LOA- We find the same qualitative features in 
all cases presented. Indeed, the discussion as in Section 
nil C[ but for sapphire, reveals that it is only for magnetic 
transitions we obtain a behavior Fg{z, uja) oc l/(fc^z)^ or 
l/{kAz)^ for sufficiently small fc^z as presented in FiglT] 
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APPENDIX A: GREEN'S TENSOR FOR A SLAB 

The equal position scattering Green's tensor for the ge- 
ometry as shown in Fig. [l] i.e. a semi-infinite half space, 
is (z > and lo is real) 

X {r,r, uj) xV = 



CnW 



+ Ca./(A) — 



??o(w) 

1 
1 
1 














772 (c.) 
2A2 



(Al) 



This equation is valid for a finite thickness h of the slab, 
i.e. it is more general than the analogous Eq. (Gl) in Ref. 
[29( 1 . Note that the analytical continuation of Eq. (jAl[) is 
Eg. ([25]) . The scattering coefficients Cat (A) and Ca/(A) 
are given by [36j 



1 



„i2r,(A)/i 



C^(A)-rp(A) Y„^2^x)e^2,iX)h 

2 _ gi2jj(A) h 



Cm(A) =r,(A) 



1 -r2(A)e'2';(A)'i ' 



(A2) 



(A3) 



with the electromagnetic field polarization dependent 
Fresnel coefficients 



'■.(A) 



(A4) 



Here we have defined 



e{uj)f]o{X) - T]{X) 
eH?7o(A)+ry(A) 



77(A) - y/kMco)-X^., 



and 



Vo 



(A) = \/k^ - A2 , 



(A5) 



(A6) 



(A7) 



as well as fc = uj/c. 

For the case of electric dipole transitions, the equal 
position scattering Green's tensor G'^(r, r, oj) rather than 
^ X G'^(r, r,a;) x ^ is needed. This Green's tensor can 
be obtained from Eq. (jAl[) by the substitution Cat (A) <-> 
Ci\i{X) and multiplying with (? juP' . In addition there will 
be a up' factor in the calculation of a transition rate F due 
to the presence of the electric field in the electric dipole 
Hamiltonian. 



APPENDIX B: PERFECT CONDUCTING SLAB 

For an excited atom above a perfect conducting slab 
as shown in Fig. [1] the frequency shift is 



(Bl) 



47rn (2z)3 

where we have defined the dimensionless functions 

/||(x) = - (2a;)2 (cos(2x) + i) 

+ 2a;sin(2a;) -f cos(2a;) — 1 

+ (2x + l)e-2- , (B2) 



and 



fi_ {x) = 2 /|| {x) ~ 2 (2x)2 cos( 2x ) . (B3) 



Here ip{kpaz) (using p = j|, _L) is given by Eqs. ([59]) and 

601). 

For short distances, i.e. kfjaZ <C 1 (non-retarded 
limit), f±{kpaz) ~ 2 f\\{kfjaz) « 2 in which case Eq. ljBip 
is reduced to the same result as for the state i.e. 
Eq. (|¥Tl) . Since for a two- level system — I/if"! this 

agrees with the intuitive fact that the CP force for the ex- 
cited (I a)) and the ground (|/3)) state should be the same 
in the classical limit, i.e. for ki^aZ <C 1. In Ref.jl^ Bar- 
ton has realized the same situation for the energy shift 
in the two-level approximation for electric dipole tran- 
sitions, i.e. the corresponding attractive short distance 
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CP force is the same in the ground state and the excited 
state. 

In the long distance retarded hmit, i.e. k/^aZ 3> 1, we 
reahze that 

/|| {k(3az) ~ - (2 kfiazf ^cos(2 kf^az) + ^ , (B4) 

and 

/j. (kp^z) ~ - {2ki3^zf (4 cos(2 kf3az) + l], (B5) 
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